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We discuss the influence of closed Dirac strings on the photon propagator in the 
Landau gauge emerging from a study of the compact £7(1) gauge model in 2 + 1 
q | dimensions. This gauge also minimizes the total length of the Dirac strings. Closed 

Dirac strings are stable against local gauge-fixing algorithms only due to the torus 
i ■ boundary conditions of the lattice. We demonstrate that these left-over Dirac strings 

CN . are responsible for the previously observed unphysical behavior of the propagator of 

space-like photons (Dt) in the deconfinement (high temperature) phase. We show 
| how one can monitor the number A3 of thermal Dirac strings which allows to separate 

the propagator measurements into A3 sectors. The propagator in A3 7^ sectors 
■ is characterized by a non-zero mass and an anomalous dimension similarly to the 

confinement phase. Both mass squared and anomalous dimension are found to be 
^ I proportional to A 3 . Consequently, in the A 3 = sector the unphysical behavior 

of the Dt photon propagator is cured and the deviation from the free massless 
propagator disappears. 

Q_; PACS numbers: 11.15.Ha, ll.10.Wx, 12.38. Gc 
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> ! I. INTRODUCTION 

•i-H . 

X' 

The interest in the three-dimensional compact electrodynamics (CQED3) has two roots: 
(i) this model shares similar features with QCD such as confinement [l| and chiral symmetry 
breaking [2| and (ii) cQED 3 has applications to condensed matter systems such as Josephson 
junction arrays and high-T c superconductors 0]. All non-perturbative features of cQED 3 
arise thanks to the compactness of the gauge field, which, in turn, leads to the appearance of 
monopoles. The monopole plasma at low temperature phase gives rise to the confinement of 
electric charges |jj whereas at high temperature the confinement disappears due to binding 
of monopoles and antimonopoles into dipoles 0, . 

The confinement property manifests itself also in the gauge boson propagator in the 
Landau gauge 0]. The effect is twofold: first, an "anomalous dimension" a appears which 
modifies the momentum dependence of the propagator, and second, a mass m is generated 
which can be well understood in terms of Polyakov's theory 1]. As it is shown in Ref. Q, 
all nontrivial effects reside exclusively in the singular fields of the monopoles. At the critical 
temperature - where the monopole plasma turns into the dipole plasma - both effects 
disappear. 

The monopole binding is observed in the zero-temperature model, too, in the presence of 
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matter fields. At sufficiently strong coupling between gauge and matter fields the monopole 
plasma also turns into a dipole plasma (a, |9( . At weak coupling the dynamical matter fields 
have some influence on the anomalous dimension of the gauge boson propagator |9l llOj]. 

Once monopoles and antimonopoles have turned into pairs they cannot contribute to non- 
perturbative effects like anomalous dimension and mass generation. Therefore it is natural to 
expect that in the high temperature phase of CQED3 the masses and anomalous dimensions 
characterizing the photon propagator have to vanish. However, the numerical results of 



Ref. 11 seem to indicate the existence of a non-zero mass and anomalous dimension for 



the propagator of the spatial photons even in the high temperature (deconfinement) phase. 



It was suggested in Ref. |11| that the spatial photons are affected by a severe Gribov copy 
problem which might lead to unphysical results in the Landau gauge. 

The fact that the U(l) gauge theory in four space-time dimensions has a gauge fixing 
problem related to the Dirac strings was discussed in Ref. [HJ . It was pointed out there that 
those gauge copies which possess so-called double Dirac sheets (DDS) give rise to a wrong 
behavior of the gauge boson propagator. 

The DDS is a classical solution in the U(l) gauge model [lij that can be considered 
as a world surface of a tightly bound pair of oppositely "charged" Dirac strings. The 
configuration containing a DDS is a Gribov copy of another configuration with zero number 
of double sheets. The DDS's are closed by wrapping around the lattice torus and are not 
associated with any monopoles. It was shown in Ref. [13( that the practical removal of the 
DDS's is a quite delicate problem. 

In the present paper we study the gauge boson propagator in high-temperature cQED 3 . 
We focus on gauge configurations containing Dirac strings looping along the shortest com- 
pactified (i.e. temperature) direction. Below we shall call these loops "thermal Dirac loops" 
(TDL). We will show that configurations without such closed Dirac loops provide physically 
sane results (vanishing anomalous dimension and mass: a = 0, m = 0) for the propagator 
in the Landau gauge whereas the presence of even a single closed Dirac loop leads to a non- 
vanishing ac and m. We do not search for the double Dirac loops of opposite orientation (a 
three-dimensional analog of the double Dirac sheet discussed in Ref. |l3j). In our approach 
those double Dirac loops would be treated as two Dirac strings. 

In Section |n] we recall the lattice model and remind the tensorial structure of the photon 
propagator at T 7^ 0. In Section imi we present the numerical results establishing the relation 
between parameters of the photon propagator and closed Dirac loops. Our conclusions are 
formulated in the last Section. 



II. THE MODEL AND ITS PHOTON PROPAGATOR AT FINITE T 

In this Section we briefly describe the model, the structure of the propagator and the 
algorithms which were used in our work. All these technical details are the same as described 
in our earlier paper An interested reader may consult that paper for a more detailed 
description. 

We use the standard Wilson single-plaquette action of cQED 3 , 

S[6 l ]=pJ2(l-cos6 p ), (1) 

V 

where Q p is the U(l) field strength tensor represented by the plaquette curl of the compact 
link field Q\. The lattice coupling constant f3 is related to the lattice spacing a and the 
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continuum coupling constant - which has dimension mass 1 / 2 - of the 3D theory as follows: 

(5 = l/{agl). (2) 

The lattice corresponding to finite temperature is asymmetric, L 2 S x L t , L t <C L s . In the 
limit L s — > oo, the temporal extension of the lattice is related to the physical temperature, 
L t = l/(Ta). Using (J2J) the temperature can be expressed in units of g 2 in the following 

way: 

T (3 



9l L 



(3) 



Thus the temperature T is proportional to the lattice coupling (3: the low-temperature 
(confinement) phase is realized at small values of f3, the high-temperature (deconfinement) 
phase corresponds to large (3. 

All our simulations are performed on a 32 2 x 8 lattice. For this lattice the phase transition 
happens at [5J| f3 c = 2.30(2). In the confinement phase the density of the monopoles is 
relatively high. The gauge dependent Dirac string bits through Dirac plaquettes defined 
below (and needed to construct the gauge independent monopoles) form either connected 
clusters of open Dirac strings with monopoles and antimonopoles at their ends or clusters of 
closed strings. Therefore, for a high monopole density the number of Dirac strings is large, 
too. 

In this paper we are interested in the quantitative effects of the temporal Dirac strings. 
Therefore we would prefer to work at low monopole densities in order to be able to easily 
separate closed and open Dirac strings unambiguously. We are going to present a quanti- 
tative demonstration that thermal Dirac loops create an unphysical behavior of the trans- 
verse photon propagator due to our attempt to fix the configuration to Landau gauge. 
Thus, we have chosen to illustrate this by a simulation at (3 — 2.6 which is located al- 
ready deep in the deconfinement phase. At this value of (3 the density of the monopoles is 
quite low, p mon = 1.5(2) ■ 1CT 4 , compared to the confinement phase (where, for example, 
Pmon = 0.1950(1) at (3 — 1.0). We have mentioned already that the construction of the Dirac 
string is gauge dependent. An ideal Landau gauge fixing makes the open strings between 
the remaining monopoles and antimonopoles (in the form of dipoles) straightened, whereas 
the closed ones, if not "wrapping" , will collapse. 

The final discussion of the photon propagator will be given in lattice momentum space. 
With the vector potential Ag jfl defined in a specified gauge, the propagator is written in 
terms of the Fourier transformed gauge potential, 

which is a sum over a certain discrete set of points x = n + |/2 forming the support of Ag tfi 
on the lattice. These are the midpoints of the links in /i direction. Here n denotes the 
lattice sites (nodes) with integer Cartesian coordinates. The propagator is the gauge-fixed 
ensemble average of the following bilinear in A, 

dap) = (Aj-O ■ 

We use the sine-definition of the gauge potential: 

A n+U,» = sin {6ft )IM ) I (g 3 a) = (U HtfM - C/i J /(2 i g 3 a) . (6) 
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The lattice momenta p on the left hand side of (J3J) are related to the integer valued Fourier 
momenta k as follows: 

pM = ~*™^, ^ = 0,±1,...,±^. (7) 

In the finite temperature case, the propagator can be parameterized by three scalar 
functions (formfactors) , 

D,M = P»MDt(\t>Ips) + P^P)D L (\plp 3 ) + ^ F{lP ^ P3) , (8) 

where we have identified the temperature direction with the \x = 3 axis, and have introduced 
the notations p 2 = p\+p\ and |p| = a/p 2 ". If the Landau gauge is exactly fulfilled, one 
would expect that F(p 2 ) = 0. In our simulations presented below this function is indeed 
very close to zero. 

Eq. (jHJ) contains also two projection operators, namely, the transverse (with respect to 
the temperature direction) projection operator P T and the two-dimensional longitudinal 
projection operator P L , respectively (with i,j = 1,2), 

f%® = k - p -f , pUp) = PIM = p?M = o , (9) 
p l ,M = pAp) - p»M , p^p) = ^-^. (io) 

In the static limit, p 3 = 0, the scalar function Dl is equivalent to the correlator of tem- 
poral photons Z?x,(|p|,p 3 = 0) = D 33 (|p|,p 3 = 0). The properties of this type of propagator 
have been discussed in Ref. |7]. Analogously, the scalar function describes the behavior 
of the spatial photons. The behavior of this formfactor Dt is of special interest in the present 
paper. 

In order to analyze the effect of the Dirac strings on the propagator in general we separate 
singular (monopole) and regular (photon) contributions to the lattice gauge field on the level 
of the link angles d x following Refs. M, Q , 

6 = 0P hot + mono , 6 mono = 2nA- 1 5p[j] , (11) 

where the dual zero-form *j represents the monopoles on the dual lattice sites (i.e. the 
monopoles are defined on the cubes of the original lattice), A -1 is the inverse lattice Lapla- 
cian. The one-form on the dual lattice p[j] defines the Dirac strings that connect monopoles 
and antimonopoles because of the condition S*p[j] = *j. 

The photon part phot is free of singularities whereas the monopole part Q mono contains 
the information about the monopole and Dirac string singularities: 

i-d[d^ hot ] 2w = , ±-d[d6~} 2n = j . (12) 

Here the DeGrand-Toussaint definition of the monopole |lf| was used. Thus, besides the 
total propagator (j3J) below we will also study the singular contribution to the propagator, 
D mono (p) = (i~^mono\ and regular contribution D^Hp) = (A phot A ph ° t ). The total con- 

v ' 1 k,fi —k,v ' x ' ' k,/x —k,v' 

tains also the mixed contribution D™ xed (p) = (A^ ot A" 2 ™ ) , which is not explicitly studied 
in this paper. 
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It turns out that the momentum dependence of the formfactors Dl and D? for ming the 
total propagator can accurately be described in both phases by the functional form 0, lid lll | , 



D(p 2 



m 



2a 



P p2(l+a) _|_ m 2(l+a) 



c. 



(13) 



where a is an anomalous dimension and m a mass parameter; Z represents the renormaliza- 
tion of the photon wave-function and C corresponds to a point-like interaction between the 
photons which we relate to a lattice artifact. In the following we denote the fit parameter 
for the D L / T formfactors by ollit etc. 

In Figure Q wepresent the fitting results for the Dl and Dt formfactors, which were 
obtained in Refs. (3,bj|, as functions of (3 (temperature). The formfactors were studied in 
the static limit, P3 = 0, and then fitted by the function ()13j) where p 2 was identified with 
p 2 . One can see from Figure ^a) that the anomalous dimension for the temporal photons, 



0.3 



a 



0.2 



0.1 - 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 



i 



_i i i i 



1.5 



D, 



* ■ i 



_i i i i i i i i i i 



2.5 



P 



(a) 




(b) 



FIG. 1: Fit parameters cxl/t and m L/T vs - P f° r At, and Dt formfactors (from Refs. |7L llll 

aL, vanishes exactly at the point of the phase transition (j3 c = 2.30(2) according to Ref. 0) 
whereas the spatial photons do not seem to feel this transition. This is documented by «t 
not vanishing at (3 = f3 c . Similarly, a different behavior is seen for the mass parameters mi 
and rriT as shown in Figure Wlb). 

The vanishing of the parameters and of Dl at (3 C and beyond is clearly corrob- 
orating the finite temperature phase transition 0] caused by dipole formation. The fields 
of the dipoles are weak at large distances and they cannot cause neither the Debye screen- 
ing [17] nor confinement 1 . Thus the origin of the non-zero mass of the spatial photon 
propagator in the deconfinement phase is likely to be an artifact of the lattice simulations, 
more specifically, of the insufficient gauge fixing. We have checked in Ref. 0] that this 



propagator is strongly affected by the Gribov copy problem. As we have argued in Ref. [11 
the Landau gauge minimizes the total length of the Dirac strings and, in the idealized case, 
does not allow for closed Dirac strings (Dirac loops) to exist. However, if the Dirac loop 



The dipole gas modifies only a short-distance interaction between electric charges providing a small linear 
correction to the Coulomb interaction [l8| . 
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is closed by wrapping around the torus, it is practically impossible to get rid of this gauge 
artifact relying only on local gauge- fixing algorithms. Gauge copies with unremoved TDL's 
correspond to local minima of the gauge fixing functional. We have qualitatively noticed 
in Ref. l3| that these strings should be blamed for the unphysical behavior of Dt in the 
deconfined phase if compared to Dl- 



III. QUANTIFYING THE EFFECT OF WRAPPING DIRAC STRINGS 

Let us now discuss the effect of closed Dirac loops on the spatial photon propagator in 
a more quantitative way. As we have already mentioned, we have simulated the model at 
P = 2.6 using a 32 2 x 8 lattice. As in our previous work we have considered two possible 
update algorithms: one purely local one (five-hit Metropolis update alternating with a 
microcanonical sweep) and another which included random offers of changing the flux in one 
of the three directions by one unit, augmented by a Metropolis acceptance check. Finally, 
in order to reduce the influence of choosing between different Landau gauge-fixed copies we 
have performed the gauge fixing procedure 100 times on random gauge copies of the same 
Monte Carlo configuration (A CO p y = 100). The evaluation of the propagator was done on the 
"best" configuration corresponding to the relative maximum of the gauge functional among 
all 101 gauge fixed configurations. 



In our previous work [ll| the fit parameters which should describe the formfactor 
have been obtained for a number A copy which was gradually increased, and even 100 Gribov 
copies were found to be insufficient for convergence. This seems to exclude the possibility 
to improve the result further by local gauge-fixing algorithms exclusively. The parameters 
for the longitudinal formfactor Dl were found to converge already for A CO p y pa 10. This 
explains why, for the present purpose, we kept iV CO p y = 100. Due to the asymmetry of the 
lattice it is relatively easy that Dirac strings are generated running around the lattice in 
temporal direction. We restrict ourselves to a string search in that direction in order to 
separate our ensemble of (locally) gauge-fixed configurations into classes according to the 
number of TDL's. A Dirac loop is formed by a connected sequence of dual links carrying 
directed Dirac string bits nS"™ ^ 0. A full Dirac string is defined by penetrating a stack 
of Dirac plaquettes. A plaquette, say P = Pg- >fW = P( Xl ,x 2 ,x 3 y,i,2, is identified as one of the 
pierced, i.e. Dirac plaquettes if 

= ^P]m0d2. • (14) 

is a non-vanishing integer. Positive or negative values - usually plus or minus one - define 
the direction of the string bit. 

The implementation of such a search, following the stack of Dirac plaquettes and counting 
the number of closed Dirac loops, is too time consuming in general. Deeper in the deconfined 
phase, it is possible to characterize a given thermalized and gauge fixed configuration by 
summing up the modules A 3 and the values of those integer Dirac plaquettes I3 pointing 
into the third ("short") compactified direction 



A3 = Integer 



-y 



, Dirac I 



x.3 



(15) 
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and 

(16) 

The different sectors of possible thermal Dirac loops are labelled by a number N% = 0, 1, . . . 
If the expected TDL's are mainly static (i.e. already minimized in length), the quantity 
N$ counts the number of those strings under the assumption that the number of monopoles 
is already very low (what is the case at (3 — 2.6). In case I/3I coincides with N$ > 1, all 
TDL's have the same wrapping orientation, otherwise TDL's with different orientation are 
present. With this procedure we cannot assess whether those strings are really "at rest" on 
the same (xi,x 2 ) position or not, but this classification is robust enough to allow for some 
local dislocations of the TDL's. The ratio 

Number of measurements with 1 1 3 1 = N 3 
All measurements in sector N$ 

counts the fraction of such strings having the same wrapping orientation. Having identified 
the "thermal Dirac loop number" N 3 , we can classify a given gauge-fixed configuration 
according to its string content and measure the photon propagator separately in the different 
sectors. 

In our previous Monte Carlo analysis |ll| we have used both local and "blended" up- 
dates, the latter including in addition to the local update adding and subtraction of fluxes, 
among them fluxes through the (x%, a^-plane. Therefore, in the light of the string sector 
classification, we also have to consider the effect of the chosen Monte Carlo algorithm on 
this classification. 

First, we have checked the dependence of the propagators on the number of thermal Dirac 
loops. In Figure 121 we present the formfactors Dl, Dt and F in the sectors without and with 
one TDL. Here we have used only the local update algorithm. One may readily notice that 
the total formfactor of the spatial photons Dt shown in Figure |2Ja,b) depends significantly 
on the number of thermal Dirac loops whereas the Dl and F formfactors are insensitive 
(within error bars) to that number. 

A very similar effect is observed for the singular contributions to formfactors which are 
depicted in Figures|2tc,d). The singular contribution to the Dt formfactor in the sector with 
one TDL (Figure Efd)) is about two orders of magnitude larger than in the N 3 = sector, 
shown in Figure I2jc). The singular contributions to other formfactors are not affected by 
the presence of TDL's. 

In order to show that the thermal Dirac loops make contributions only to the singular 
and mixed part of the Dt propagator, whereas the regular part is not affected by the Dirac 
string, we plot in Figures Efe,f) the regular part of the Dt, Dl and F formfactors as a 
function of spatial momentum. One can see that the regular contributions in the zero-loop 
and one-loop sectors coincide within error bars. 

Qualitatively, the results in Figures El can be understood as follows. The fact that the 
regular contribution to the propagator is insensitive to the number of TDL's is very natural 
since the regular part does not receive contributions from singular structures like the Dirac 
strings. The sensitivity of the Dt formfactor to the number of TDL's and the respective 
insensitivity of the Dl formfactor follows from Eq. (JHJ. Indeed, the TDL is described by a 
chain of the plaquettes which are perpendicular to the temporal direction. The boundary of 
the plaquette 5p[j] is a vector field the spatial components of which are non-zero, whereas 
the temporal ones are vanishing. The inverse Laplacian in Eq. (fTTj) does not mix these 
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FIG. 2: The formfactors Dl, Dt and F vs. spatial momentum in the sector with (left column) 
zero, N3 = 0, and (right column) one, N3 = 1, thermal Dirac loop using the local update. The 
total (a,b), singular (c,d) and regular (e,f) contributions are shown. 
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components. Thus, a thermal Dirac loop may provide a contribution only to the spatial 
components of the gauge potential, i. e. to only. Finally, the insensitivity of the formfac- 
tor F to the presence of TDL's follows from the fact that F corresponds to the longitudinal 
(in momentum) part of the propagator whereas the contribution from the TDL is transverse 
(see Eq. (JTTJ)). 

In the sector without thermal Dirac loops (N3 = 0) the tiny singular contribution to the 
propagators Dt and Dl can be explained as an effect of remaining monopole-antimonopole 
pairs which are mainly oriented in temporal direction (and/or additionally of strings in 
spatial direction on which we did not trigger). 

So far, the results in Figures 121 were obtained from the evaluation of gauge-fixed config- 
urations when the original ensemble was generated by exclusively local updates. We have 
repeated the same analysis for the blended update algorithm mentioned which includes also 
global changes of fluxes. On this basis, we notice that (within error bars) the results for the 
propagator formfactors in the individual TDL sectors are independent of whether blended 
updates are allowed or not. 

However, the choice of the update algorithm (blended or local) influences the relative 
weight of the individual sectors within the gauge-fixed ensemble. Strictly speaking, the 
content of wrapping Dirac strings is a gauge artifact. We are only monitoring it. Only 
dedicated, "big" gauge transformations would be able to remove them. For example, if 
global flux changes are accepted in the blended update then after applying local gauge fixing 
the sector with N 3 = ±1 is slightly dominating. If a local update is applied without offering 
global flux changes, after local gauge fixing the N 3 = sector is clearly the dominating one. 
An overview of the statistics available for this study is given in Table [Q 



iV 3 (Update) 


r3 


# of meas. in sector ./V3 


# of meas. attempts 


(local) 


1.0 


1191 


2000 


(blended) 


1.0 


767 


2000 


1 (local) 


0.891 


911 


2500 


1 (blended) 


0.903 


1115 


2500 


2 (local) 


0.866 


1912 


7000 


2 (blended) 


0.778 


974 


7000 


3 (local) 


0.859 


580 


15000 


3 (blended) 


0.869 


314 


20500 


4 (local) 


0.895 


19 


600000 


4 (blended) 


0.931 


72 


600000 



TABLE I: Statistics for the different thermal Dirac loop sectors using local or blended update. 

From this Table we can speculate that in the process of gauge fixing the local algorithm 
might have annihilated thermal Dirac loops of different orientation to a large extent. This 
is expressed by the ratio (fT7|) which is of the order of 80 % or larger. This would correspond 
to a suppression of so called double Dirac loops (analogues of the mentioned DDS in four 
dimensions) in the configurations possessing two oppositely oriented TDL's. 
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It turns our that the fit of the formfactor using the functional dependence on p given 
in f|13|) works very well separately for all sectors iV 3 = 0, ... ,4, with x 2 /d.o.f. ~ 0.5. The 
corresponding fitting curves yielding a? and tut separately for each jV 3 = 0, . . . , 3, are shown 
in Figure El 
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FIG. 3: The total formfactor Dt for spatial photons in the different TDL sectors, using the blended 
update. The data are fitted by the function (|T3)). 

The fitting parameters themselves as functions of the TDL multiplicity iV 3 are shown in 
Figure |U both for the local and the blended updates. First we notice that the fit parameters 
in a given loop sector are practically independent of the update. As one can see from Table HJ 
our statistics for iV 3 = 4 achieved by exclusively local updates is quite low such that we omit 
the corresponding points in Figures 0] We emphasize the observation that it is extremely 
unlikely to produce, by updates without global changes of flux, configurations which finally, 
after gauge fixing, end up in the sector iV 3 = 4. 

In the sector without thermal Dirac loops (N 3 = 0) the anomalous dimension cut, the 
mass parameter ttit and the contact term parameter Ct are consistent with zero. At the 
same time, the renormalization of the photon wavefunction Zt is very close to unity. Thus, 
if we would restrict ourselves to the iV 3 = sector in the propagator measurements, there 
is practically no difference between the Dl and Dt formfactors in the deconfinement phase, 
and we reproduce the expected behavior of a trivial free photon propagator. 

In the non-zero TDL sector (iV 3 ^ 1) the parameters a?, and Ct become non- 
zero, and the parameter Zt deviates from unity. Thus, the formfactor acquires a non- 
trivial momentum dependence compared to the sector without thermal Dirac loops. The 
dependence of the propagator parameters, «t, m T , Z T and Ct, on the multiplicity of thermal 
Dirac loops, iV 3 , can be described by simple linear functions: 



olt = da A3 



rrin 



,.A, 



1 = a z N 3 , C T = a c N 3 



(18) 



The corresponding fits are shown by the dashed lines in Figures EJ The quality of these fits is 
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FIG. 4: Fit parameters ay, tut, Zt and Ct for the Dt formfactor vs. the number of thermal Dirac 
loops N$ using data obtained with blended or local updates. The fits by Eqs. (|18|) are presented 
by dashed lines. 

approximately the same, x 2 /d.o.f. ~ 1.5. The proportionality coefficients are a a = 0.062(5), 
a m = 0.0745(38), a z = 0.84(3) and a c = -0.011(1). We recall that all this refers to a 
temperature corresponding to (3 = 2.6. We expect that in general these coefficients must 
depend on temperature. 



IV. CONCLUSIONS 

We have studied the influence of thermal Dirac loops on the gauge boson propagator 
within the finite temperature compact U(l) gauge model in three dimensions. We have used 
the Landau gauge and worked deep inside the deconfinement phase where monopoles are 
very dilute and Dirac strings are dilute as well. This allows us to unambiguously recognize 
the thermal Dirac loops. Wrapping Dirac strings are ubiquitous on a finite lattice. On an 
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asymmetric lattice at high temperature they are predominantly closed along the temperature 
direction which is the shortest lattice direction. Although closed Dirac loops along the spatial 
directions (non-thermal Dirac loops) are not completely excluded, they are extremely rare 
and do not yield an essential contribution to the propagator D T in question. 

Strictly speaking, in the Landau gauge Dirac strings closed due to periodic boundary 
conditions are artifacts of the gauge fixing, because the Landau gauge condition corresponds 
to the minimization of the number and length of Dirac strings We have found that the 
presence of such thermal Dirac loops seriously affects the properties of the propagator, in 
the considered case those of D?, in the deconfinement phase. The propagator formfactor 
corresponding to spatial photons in a sector with a non- vanishing number A 3 of TDL's 
mimics a momentum dependence similar to what is known from the confinement phase, 
Eq. (|13p. The parameters, which describe the deviation from the free photon propagator, 
which ought to be expected in the deconfinement phase, are found clearly proportional to 
N 3 because we were working in dilute gas regime. An explanation of this coincidence may 
lie in the fact that the monopoles, which are active in the confinement phase, contribute to 
the propagator only indirectly, i. e. via the Dirac strings. 
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